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Abstract 



This paper is about the apphcation of geometric algebra to black 
holes. Geometric algebra is the natural extension of the concept of 
multiplication from real numbers to geometric objects. It is a Clif- 
ford algebra with a geometric interpretation. The most important 
tool within geometric algebra used here is the Gauge Theory of Grav- 
ity (GTG). The theory uses gauge fields to ensure that all relations 
between physical quantities are independent of the position and orien- 
tation of the matter fields. These are position-, rotation-, and phase- 
gauge fields. 

Applications begin with a preliminary discussion on cosmological 
Hawking radiation and the Unruh effect. The result for the 'Hawking 
temperature' of the de Sitter cosmological horizon is quite interesting. 
If a thermodynamic system has a negative temperature then that sys- 
tem should be unstable, and tend to expand or collapse. But, if we in- 
stead abandon the concept of particles as being observer-independent 
in this case, we find the usual positive 'Hawking temperature'. Be- 
sides, geometric algebra was also successfully applied to treat the Un- 
ruh effect. 

The discussion then moves on to consider rclativistic wave equa- 
tions in spherically symmetric static black hole backgrounds. The 
Klein-Gordon and Dirac equations are treated on Schwarzschild and 
Reissner-Nordstr5m black hole metrics. The analysis is generalised 
to allow for the presence of magnetic monopoles. An interesting sit- 
uation concerning electromagnetism in the Schwarzschild black hole 
background is treated. Then it is repeated on the axially symmetric 
stationary (i.e., Kerr) black hole background. The Hawking temper- 
atures were derived for all cases and compared with the results of 
'standard' calculations. Remarkably, they agreed. 

It was also found that the absolute values of the real part of the 
non-zero roots in the Klein-Gordon, Dirac, and electromagnetic field 
cases are the values of the spin of the corresponding particles yet 
electromagnetism, at least, is completely classical. It seems that geo- 
metric algebra can tell that Maxwell's purely classical theory concerns 
spin-1 particles (photons) even though the picture of bosonic particles 
carrying forces is a quantum description! In some sense, the quantal 
'geometric' nature of the theory is recognised by geometric algebra! 



1 Introduction 



The study of black holes began with the classic paper of Oppenheimer and 
Snyder [221 • They considered the collapse of a homogeneous spherical gas 
without pressure in Einstein's General Relativity (GR). They found that the 
sphere eventually becomes cut off from all communication with the rest of 
the universe, and at the end the matter was crushed to infinite density at the 
centre. Black holes were largely ignored until Wheeler and his collaborators 
revived the subject in the 1960's and triggered the golden era of black hole 
research. 

GR clearly predicts the existence of black holes, which are the most ex- 
citing predicted compact objects in the universe (see, for example, Kramer, 
Stephani, MacCallum, and Herlt |40|, Misner, Thorne, and Wheeler 
Hawking and Ellis |31J, Christodoulou [TU], and de Felice and Clarke [T^ : 
for popular accounts see, for example, Einstein Thorne |Snj, Begelman 
and Rees jl]. Wheeler and Ford jHE] and Gribbin [2H])- The Gauge Theory 
of Gravity (GTG) employed in this dissertation also predicts black holes as 
in GR but with some modifications as found in Lasenby, Doran, and Gull 
|41j . These solutions have spacetime singularities from which nothing can 
emerge to infinity. There is also a null hypersurface called the event horizon 
or 'surface' of the black hole, which separates the spacetime within it from 
the rest of the universe. In Newtonian language, the escape velocity from 
inside the surface exceeds the speed of light, so that nothing can escape. The 
concept of the event horizon is purely classical. If the singularities are hid- 
den behind the event horizon, nothing goes wrong with GR or GTG, because 
anything which comes in cannot go out, so we can argue that the ingoing 
information is stored behind the event horizon even if we cannot access it. 
But, if Quantum Mechanics (QM) is involved, it seems that GR or GTG 
will break down because smooth initial data (pure states) are able to evolve 
into singular field configurations (mixed states). Hawking showed semiclas- 
sically that QM implies the emission of particles from black holes, which is 
classically impossible [3U] (see also and [M]; see |^ for the interpretation 
of Hawking radiation as tunneling). The radiation which escapes is purely 
thermal, as if a black hole is a black body with a temperature proportional 
to its surface gravity and carries no information at all. When a black hole 
evaporates completely, all the information of the collapsed matter is lost. 
This phenomenon, which is usually called information loss or quantum de- 
coherence, would imply non-unitary evolution in contradiction to one of the 
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basic laws of QM. Instead of discussing the information loss phenomenon, we 
shall consider below how the Hawking radiation appears naturally in GTG 
within the language of geometric algebra. 

This paper is mainly about the applications of geometric algebra to black 
holes. Geometric algebra is firstly introduced. This powerful mathematical 
tool extends the concept of multiplication from real numbers to geometric 
objects. It is a Clifford algebra with a geometric interpretation. Its power 
to unify different languages in mathematical physics has been rediscovered 
time after time. It promises one common mathematical language through- 
out physics. In the Dirac equation, for instance, we can replace the usual 
imaginary scalar i with a quantity interpretable as rotation in a particular 
plane. This permits the disentanglement of geometry and quantum effects, 
clarifying the latter. 

The language of geometric algebra is used to express gravity as a gauge 
theory (Gauge Theory of Gravity, or GTG). The theory uses gauge fields to 
ensure that all relations between physical quantities are independent of the 
position and orientation of the matter fields. The discussion of GTG begins 
with the gauge fields. These are position-, rotation-, and phase-gauge fields. 

An interesting example on the cosmological background, i.e., the de Sitter 
metric, that relates to cosmological Hawking radiation and the Unruh effect 
is given. 

The relativistic wave equations on the spherically symmetric static black 
hole backgrounds are then discussed. The Klein-Gordon and Dirac equations 
are treated on the Schwarzschild and Reissner- Nordstrom black hole metrics. 
The Hawking temperature for each case is derived and compared with the 
result from the 'standard' calculations. Generalisation to allow for the pres- 
ence of magnetic monopoles is also presented. An interesting case concerning 
electromagnetism in the Schwarzschild black hole background is examined as 
well. 

The relativistic wave equations on the axially symmetric stationary black 
hole backgrounds are investigated. The Klein-Gordon and Dirac equations 
are treated on the Kerr black hole metric. The Hawking temperature for 
each case is again derived and compared with the result from the 'standard' 
calculations. 

Finally, there are two appendices which comprise short notes on spherical 
monogenic functions and vector spherical harmonics. The computations in 
this paper is performed with the help of Maple package |57j . 
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2 Geometric Algebra 



This brief introduction to geometric algebra is intended to establish out no- 
tation and conventions. More complete introduction may be found in |1T] 
and [32j. The basic idea is to extend the algebra of scalars to that of vectors. 
We do this by introducing an associate (Clifford) product over a graded lin- 
ear space. We identify scalars with the grade elements of this space, and 
vectors with the grade 1 elements. Under this product scalars commute with 
all elements, and vectors square to give scalars. If a and b are two vectors, 
then we write the geometric (Clifford) product as the juxtposition ab. This 
product decomposes into a symmetric and an antisymmetric part, which de- 
fine the inner and outer products between vectors, denoted by a dot and a 
wedge, respectively : 

a.b = ^{ab + ba) 
a Ab = |(a6 - ba) ^ ^ 

It is simple to show that a.b is a scalar, but aAb is neither a scalar or a vector. 
It defines a new geometric element called a bivector (grade 2). This may be 
regarded as a directed plane segment, which specifies the plane containing a 
and b. Note that if a and b are parallel, then ab = ba, whilst ab = —ba for a 
and b perpendicular. This process may be repeated to generate higher grade 
elements, and hence a basis for the linear space. 

2.1 The Spacetime Algebra (STA) 

The Spacetime Algebra is the geometric algebra of spacetime. This is familiar 
to physicists as the algebra if the Dirac 7-matrices. The STA is generated 
by four orthogonal vectors {7;^}, /i = . . . 3, satisfying 

lf^-'lu = ^{ifilu + lul^,) = Vfiu = diag(+ ). (2) 

A full basis for the STA is given by the set 

1 {7^} Wk,i(^k} i 

1 scalar 4 vectors 6 bivectors 4 trivectors 1 pseudoscalar (3) 
grade grade 1 grade 2 grade 3 grade 4 

where = 7^70, k = 1 ... 3, and i = 7o7i7273 = cricr2cr3- The pseudoscalar 
i squares to —1 and anticommutes with all odd-grade elements. The {ck} 
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generate the geometric algebra of Euclidean 3-space, and are isomorphic to 
the Pauli matrices. They represent a frame of 'relative vectors' ('relative' to 
the timclike vector 70 used in their definition). The {(Jk} arc bivectors in 
four-dimensional spacetime, but 3-vectors in the relative 3-space orthogonal 
to 7o. 

An arbitrary real superposition of the basis elements (3) is called a 'mul- 
tivector', and these inherit the associative Clifford product of the {7^} gen- 
erators. For a grade-r multivector Ar and a grade-s multivector Bs we define 
the inner and outer products with 

Ar • Bs= {ArBs)\r-s\: A^ A Bg = {Aj-Bs) \r+s\: (4) 

where {M)r denotes the grade-r part of M. We will also use the commutator 
product 

AxB=hAB-BA). (5) 
The operation of reversion, denoted by a tilde, is defined by 

{AB)=BA (6) 

and the rule that vectors arc unchanged under reversion. We adopt the 
convention that in the absence of brackets, inner, outer and commutator 
products take precedence over Clifford products. 

Vectors are usually denoted in lower case Latin, a, or Greek for basis 
frame vectors. Introducing coordinates {x^{x)) gives rise to a (coordinate) 
frame of vectors {e'^}, satisfies • e'' = 5^. The vector derivative V(= is 
then defined by 

V = e^d^ (7) 

where = d/dx^. 

Linear functions mapping vectors to vectors are usually denoted with an 
underbar, /(a) (where a is the vector argument), with the adjoint denoted 
with an overbar, /(a). Linear functions extend to act on multivectors via 
the rule 

/(a A 6 A • • • A c) = /(a) A /(6) A • • • A /(c), (8) 



4 



which defines a grade-preserving hnear operation. In the STA, tensor objects 
are represented by Unear functions, and all manipulations can be carried out 
in a coordinate-free manner. 

All Lorentz boosts or spatial rotations are performed with rotors. These 
are even-grade elements R, satisfying RR — 1. Any element of the algebra, 
M, transforms as 

M ^ RMR. (9) 

A general rotor may be written as R — exp(S/2) where S is a bivector in 
the plane of rotation. 

2.2 Gauge Theory of Gravity (GTG) 

Physical equations, when written in the STA, always take the form 

A{x) = B{x) (10) 

where A{x) and B{x) are multivector fields, and x is the four-dimensional 
position vector in the (background) Minkowski spacetime. We demand that 
the physical content of field equations be invariant under arbitrary local 
displacements of the fields in the background spacetime, 

A{x) ^ A{x'), x' = /(x), (11) 

where f{x) a non-singular function of x. We further demand that the physical 
content of the field equations be invariant under an arbitrary local rotation 

A(x) ^ RA(x)R, (12) 

with R a non-singular rotor- valued function of x. These demends are clearly 
equivalent to requiring covariance (form-invariance under the above trans- 
formation) of the field equations. These requirements are automatically sat- 
isfied for non-derivative relations, but to ensure covariance in the presence 
of derivatives we must gauge the derivative in the background spacetime. 
The gauge fields must transform suitably under (local) displacements and 
rotations, to ensure covariance of the field equations. This leads to the in- 
troduction of two gauge fields : h{a) and fi(a). The first of these, h{a), is a 
position-dependent linear function mapping the vector argument a to vectors. 
The position dependence is usually left implicit. Its gauge-theoretic purpose 
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Table 1: Symmetry transformations of the gravitational gauge fields 



gauge transformation 


transformed fields 




h\a, x) 


Q'{a, x) 


displacements 




h{T\a),x') 


nii{a),x') 


spacetime rotations 


Rip{x) 


Rh{a, x)R 


Rn{a,x)R-2a-VRR 


phase rotations 




h{a, x) 


Q{a, x) 



is to ensure covariance of the equations under arbitrary local displacements 
of the matter fields in the background spacetime jUj. The second gauge 
field, fl{a), is a position-dependent linear function which maps the vector a 
to bivectors. Its introduction ensures covariance of the equation under local 
rotations of vector and tensor fields, at a point, in the background spacetime. 

Once this gauging has been carried out, and a suitable lagrangian for the 
matter fields and gauge fields has been constructed, we find that gravity has 
been introduced. Despite this, we are still parameterising spacetime points 
by vectors in a flat background Minkowski spacetime. The covariance of the 
field equations ensures that the particular parameterisation we choose has 
no physical significance. The feature that is particularly relevant to this is 
that we still have all the features of the fiat-space STA at our disposal. A 
particular choice of parameterisation is called a gauge. Under gauge transfor- 
mations, the physical fields and the gauge fields will change, but this does not 
alter physical predictions if we demand that such predictions be extracted 
in a gauge- invariant manner. List of symmetry transformations which make 
the action invariant is given in Table 1, while the conventions of GTG are 
shown in Table 2 |41j . 

The covariant Riemann tensor TZ{a A b) is a linear function mapping 
bivectors to bivectors. It is defined via the field strength of the Q{a) gauge 
field : 

nh-\a Ab) = a- Vn{b) - b ■ Vn{a) + fi(a) x n{b), (13) 

The Ricci tensor, Ricci scalar and Einstein tensor are formed from contrac- 
tions of the Riemann tensor |1T] : 
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Table 2: Definitions and conventions 



displacement gauge fields 


h{a) 


rotation gauge fields 


Q{a), uj{a) = Qh{a) 


spinor derivatives 


Daip = a ■ W + \^{a)'4) 

a- Dip = a- h{V)ip + |u;(a)V' 


observable derivatives 


VaM = a ■ VM + n{a) x M 

a ■ VM = a ■ h{V)M + uo{a) x M 

VM = daa ■ VM = V- M + VAM 


vector derivative 




multivector derivative 


dx = E A ■ ■ • A eJ (ej A ■ ■ ■ A Ci) * (9x 

i<-<j 



Ricci Tensor: 7^(6) = da ■ 7^(a A b) 
Ricci Scalar: n = da ■ n{a) (14) 
Einstein Tensor: Q{a) = lZ{a) — ^aTZ. 

The Einstein equation may then be written as 

g{a) = KT{a), (15) 

where T{a) is the covariant, matter stress-energy tensor. The remaining field 
equation gives the f2-function in terms of the /i-function, and the spin of the 
matter field (41j . 

3 Cosmological Hawking Radiation and the 
Unruh Effect 

3.1 Dirac equation on the de Sitter metric 

In this section we discuss the Dirac equation on the de Sitter background 
ITTj . We employ the minimally coupled Dirac equation from ITS] . 
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which is in this case 



Dtp I as = mp^jo. (16) 
We use the Kerr-Schild gauge : 

h{a) = a + A(r)a • e_e_ (17) 

where 

A(r) = ^r^ (18) 
o 

and e_ = 70 — e^. The Dirac equation now becomes 



Viplas - [(A(r)(-7o + 7r)(9t - dr)]ipI(T-i + [(-70 + 7r)(^^)]V^«o-3 = mpip'jo. 

(19) 

Upon multiplying this equation by 70 from the left, and using the symbol j 
to represent right-sided multiplication by la^, i.e., jijj = ipla^, we have 



jdti^ = -jV^ - jA(r)(l + ar.)idt - 9,)^' + j(^r)(l + a,)^ + rUp^P (20) 
where ip = 7oV'7o- By making a trial separation of variables 

^{x) = a(t)^(x) (21) 

we find that 

a{t) = exp{-jEt), (22) 

and 

V^(x) -A(r)(l + (T,,)5,,^(x) 
-(§r)(l + a,)^(x) + A(r)(l + a,)(-jE)^(x) (23) 
-jEip{x) +jmpip{x) = 0. 

Next, we need to separate out the angular dependence. This can be done 
by using the spherical monogenics [20]. Note that, from equation (24) to 
equation (27) below, we temporarily employ Icr^, instead of j; these after we 
use j for 1(73 again. The (unnormalized) monogenic is defined by 
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: [(/ + m + l)P,'"(cos^) - Pr+^(cose)/a0]e'"*^"^ (24) 

where / > 0, — (/ + 1) < m < /, and are the associated Legendre polyno- 
mials (see, for example, PP). Two important properties of the %p]^ are 



-(//r)a,C 



and 



V(a.O = (/ + 2)/rC. 
After some manipulation [20], the solutions are as follows : 

V'(x, k) = 



ar'ipruir)a3 + ^["/^(r), k = -(/ + 1), 



(25) 



(26) 



(27) 



where «; is a non-zero integer and u{r) and v{r) are complex functions of r. By 
substituting equation (27) into equation (23) and employing the properties of 
the spherical monogenics above, we obtain the coupled equations (in matrix 
form) 



B 



Ui 
U2 



where 



(28) 



and 



B 



-A(r) 
-A(r) 



-A(r) 
1 - A(r) 



f + C(r) 3{E + mp) + ar) 



rrir. 



C(r) 



'F + C(r) 



Here A(r) is given by equation (18), and 



Afr) A 

ar) = ^-1+ r + jEX{r) 

r 6 



(29) 



(30) 



(31) 



where ui and U2 are the reduced functions defined as in the Schwarzschild 
case : 
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Ml = ru, U2 = jrv (32) 

and the primes in (28) denote differentiation with respect to r. Equation 
(28) can also be written as 



A 



u 



AB 



Ui 
U2 



with 



A 



1 - A(r) A(r) 
A(r) 1 - A(r) 



(33) 



(34) 



We now look for power series with a limited radius of convergence around 
the horizon, by introducing the series 



where 



and 



k=0 k=0 



(35) 



(36) 



(37) 



As in the Schwarzschild case the index s controls the radial dependence 
of ip at the horizon, so that it represents a physical quantity. To find the 
values of s we substitute (35) into (33) and set r/ = 0. The resulting indicial 
equation is 



det 



AB + (-V3A)sJ 



0, 



(38) 



r=r„ 



from which we find the solutions for the two indicial roots : 



(39) 



Since s = is always a solution of the indicial equation, solutions that are 
analytic at the horizon always exist [H]. In j3T] it is claimed that one can 
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calculate reflection coeflicients and scattering amplitudes by determining the 
split between ingoing and outgoing states of these solutions. However, the 
problem here is whether the second root can be physically signiflcant. We 
now consider the 'Hawking temperature', which will be denoted by Th- In 
these calculations we shall need the non-covariant conserved current (which 
generates the streamlines whose curves are timelike) that can be expressed 
as gU : 



This satisfles the flatspace conservation equation V- J = 0. We express rj^ as 



J = Uijjoi^) det{h) 



1 



(40) 




(41) 



We can now write 




(42) 



with the choice of argument 




r > Th 
r < Th- 



(43) 



We can split E into real and imaginary parts as 



If we now take the limit r 



E = Er+je. (44) 

/ o 

> W from above and below we flnd that the 



7o component of J is given by 




(45) 



where i?i(^,0) is a positive-deflnite flnite term. 
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We now calculate the radial component of J 



e'-J 




(46) 



where B2{6,(j)) is a positive-definite finite term. 

From these results we can derive the 'Hawking temperature'. By taking 
the ratio of the inward fiux to the total fiux, which is 



exp 



2n^ 



-Er 



(47) 



+ 1 



we identify a Fermi-Dirac distribution with 'Hawking temperature' given by 

1 



In the calculations above we used the inward fiux, not the outward fiux, 
because we are concerned with an observer inside the de Sitter universe. We 
therefore obtain the same result as by the standard approach, but much more 
economically. 

Consider now an observer outside the de Sitter horizon (i.e., in region II 
of Figure 2 in [2Zj)- In this case we take the ratio of the outward fiux to the 
total fiux, 



1 



exp 



(49) 



+ 1 



to obtain a Fermi-Dirac distribution with a negative 'Hawking temperature' 

1 



H 



2-Kk 



B 



(50) 



This is a remarkable result. An observer inside the de Sitter horizon (i.e., 
in region I of Figure 2 in j2Zj) will detect a positive Hawking temperature, 
while an observer outside it (i.e., in region II of Figure 2 in j^ZI) will detect 
a negative Hawking temperature ! It seems that the Hawking temperature is 
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observer-dependent. In fact this is because our concept of particle is observer- 
dependent, as explained below. 

We know that the de Sitter universe is unstable and forever expanding 
|17j . The inflationary era which is believed to have taken place very early in 
the evolution of the universe has a similar nature to the de Sitter universe 
(the difference is that inflationary era took place for a very short period 
of time, not endlessly). The negative temperature above may reflect this 
instability with respect to an observer outside the de Sitter universe. If 
negative temperature is really an indication that a thermodynamical system 
is unstable, then the Cauchy horizons (i.e., inner horizons) in black holes 
should be unstable and should not therefore exist (from the point of view of an 
observer outside the black hole). Instead of expanding like the cosmological 
horizon in the de Sitter case, we find that the Cauchy horizons actually tend 
to collapse fH] . 

We also know that inside a black hole there are particle states which 
have negative energy with respect to an external stationary observer. As 
discussed by Gibbons and Hawking [27], a cosmological event horizon has 
many similarities with a black hole event horizon. Hence, if we abandon the 
concept of particles as being observer-independent and apply this view to our 
case, in which we are concerned with the observer inside the event horizon 
of a de Sitter universe, then we get the positive 'Hawking temperature' as in 
the calculations above. 



3.2 The Unruh Effect 

Suppose a detector is being uniformly accelerated in Minkowski spacetime. 
According to Unruh [HI], this detector behaves as if it were placed in a 
thermal bath of 'real' particles with temperature (i.e., Unruh temperature) 
given by 

Tu = (51) 

where a is the acceleration of the detector. The ordinary vacuum state in 
Minkowski spacetime from the viewpoint of an accelerating observer has 
thermal properties similar to those due to true particle creation by a black 
hole as seen by a stationary observer at some constant distance from the 
black hole. Hence, accelerated systems find themselves in a thermal bath. 
This is called the Unruh effect. Both Hawking and Unruh temperatures can 
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be related to information loss, associated with real and accelerated observer 
horizons respectively (see, for example, [^21 ESI IIHI IM])- 

In an interesting paper ^1], a connection was sought between surface 
gravity and the Unruh temperature, with the aim of establishing the principle 
of equivalence between constant acceleration and 'true' gravity effects by 
globally embedding curved spaces in higher dimensional flat spaces. The 
relevant acceleration in the fiat spaces gives the correct Hawking temperature 
in the curved ones. 

Consider now a constant electric field Ec produced by two oppositely 
charged plates facing each other. If the electric field is strong enough then 
the vacuum located between the plates produces real charged particles, as 
near to the event horizon of the black hole. Each particle (with mass M and 
charge e) propagates with a uniform acceleration given by 

a ^ f . (52, 

In the Unruh effect, the propagation of the detector is fully described by a 
given classical trajectory; only the internal transitions accompanied by the 
emission or absorption of quanta of the radiation fields are treated quantum- 
mechanically. The analysis is in some sense semiclassical, like Hawking's 
treatment of black hole evaporation jHO]- However, Unruh used second- 
quantised quantum field theory. In this paper only first-quantised quantum 
theory is used. 

Unruh found that an accelerated observer in Minkowski spacetime (in 
the normal vacuum state) can detect and absorb particles. According to an 
inertial observer this absorption appears to be emission from the accelerated 
observer's detector. Similarly, an observer at a constant distance from a 
black hole detects a steady flux of particles coming out from the black hole, 
with a thermal spectrum, while an observer who falls freely into the black 
hole does not see particles radiating. This is a consequence of the principle 
of equivalence. 

We now discuss the Dirac equation on the Rindler background to study 
the Unruh effect [HS] • We still employ the minimally coupled Dirac equation 
jm HHI , which is in this case 

Diplas = iripip'jo. (53) 
We use the Kerr-Schild gauge 
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h{b) ^ b + X{z)b ■ e-e^ (54) 

with e_ = 7o — 73, and X{z) = az, where a is the acceleration. With the 
ansatz 

iPix) = {u{z)+ja3v{z))e-^^\ 
the Dirac equation becomes (in matrix form) : 



(55) 



Ui 
U2 



where 



1 - A(^) -\{z) 
-X{z) 1 - X{z) 



and 



B = 



[j{E-mp) + az) -l + Ciz) ) 



Here ({z) is given by 



Ciz)^~+jEXiz) 



(56) 



(57) 



(58) 



(59) 



and ui and U2 are the reduced functions defined as in the Schwarzschild case: 

ui — zu, U2 — jzv, (60) 

and the primes in (56) denote differentiation with respect to z. Equation 
(56) can also be written as 



(61) 



where 



.'^(1-X{z) X{z) 
- I A(;.) 1 - A(^) 



(62) 



As usual, we introduce the series 
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oo 



k=0 k=0 

with 

7] = z- zh, (64) 

and 

ZH = ^, (65) 

to generate the indicial equation 



det 



A B - 2asl 



= 0, (66) 

z=zh 



The sohitions for the two indicial roots are 

. = 0,-i + ^. (67) 

As in the previous section, the fact that s = is always a solution of the 
indicial equation implies that analytic solutions always exist at z^, which in 
this case acts as an event horizon. (An interesting account of this is given 
in Chapter 4 of Birrell and Davies [Hj-) As before, we wish to investigate 
whether the second root can be physically significant. We again express 77^ 
as 



We can now write 



ln(.-l) = ln|(.-l)|+jarg(z-l), (69) 
with the choice of argument 

/ 1 \ I 0, z > Zh 

aigiz ) = <^ ' (70) 

2a \ -TT, z <zh. ^ ^ 

Before proceeding further, we split E into real and imaginary parts as 

E = Er- je. (71) 
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If we now take the limit z - 
7o component of J is given by 



1m 



from above and below we find that the 



7o-J = 5ie~2^* I 
1, 

exp 



where is a finite term. 

We now calculate the 2;-component of J : 



z > Zh, 
z < Zh, 



(72) 



73- J = 



-2et 



exp 



2'KEr 



z > Zh, 
, z < Zh, 



(73) 



where B2 is finite. 

From the above equation we are able to derive the 'Unruh temperature'. 
The ratio of the 'outward flux' to the total flux is 



73 •-'■4 



73- J+ - 73- J_ 



exp 



27r£;,. 
a 



+ 1' 



which is, remarkably, a Fermi-Dirac distribution with temperature 



Tr, 



(74) 



(75) 



as we found in the standard approach. Thus, the ordinary Minkowskian 
vacuum state is seen by an accelerating observer to have thermal properties 
very similar to the thermal effects resulting from true particle creation by a 
black hole. The temperature (in cgs units) is 



Tr/ ^ 4 X 10"^^ a K, 



(76) 



so this effect is too small to be perceived by an ordinary laboratory detector. 

In regard to the alternative explanation at the end of the previous section, 
if we naively take Tjj = Th as suggested in [H], then we have 



(77) 
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From the discussion above we may conclude that an accelerating observer 
with acceleration in Minkowski spacetime perceives a thermal bath with 

temperature /{27ikB), while an merim/ observer in the de Sitter universe 

also detects the same temperature. In the former case the acceleration causes 
the same effect as that seen by an inertial observer immersed in a thermal 
bath with corresponding temperature, while in the latter case the creation 
of particles is described as a consequence of cosmic expansion. 



4 Spherically Symmetric Black Holes 

4.1 Dirac equation on Reissner-Nordstrom black hole 
metric 

In this section we discuss the Dirac equation on the Reissner-Nordstrom 
black hole background [211 123] ; see also [21] . A similar discussion for the 
Schwarzschild background has been given by Lasenby, Doran, and Gull [H] 
and is a special case of the present analysis for which the charge q is zero. 
We use the minimally coupled Dirac equation from [HI HH] (see also 
jnni, and jSS]), which is 

Dtplas — eAip = mpip'^Q^ (78) 
and the Kerr-Schild gauge ^ (see also |23 EH 113 EDI ) 



h{a) = a + \{r)a ■ e_e_ (79) 



where 



A(0 . ^ - ^. (80) 

and e_ = 7o — e^. 

With the gravitational field of the Reissner-Nordstrom black hole back- 
ground expressed in spherical polar coordinates, the Dirac equation becomes 



M 

VV'/cr3-e^V^-[(A(r)(-7o+7r)(9t-9r)]V^/(T3+[(-7o+7r.)(^)]V'«o-3 = "^pV'7o- 

(81) 
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By multiplying this equation by 70 from the left and using the symbol j to 
represent right-sided multiplication by la^, i.e., jip = iplas, we find 



M — 

jdtijj = - j - J A(r) (1 + ar) (dt -dr+ je$)^ + i(^) (1 + ^r)ip + 

(82) 

where ip = 7oV'7o; and $ is the potential of the black hole event horizon. We 
make a trial separation of variables 

tp{x) = a(t)^(x) (83) 
which gives as the solution of the t-equation 

a{t) = exp{-jEt). (84) 
For the spatial dependence we obtain 

VV^(x) -A(r)(l + (T,)a,V^(x) 

-(^)(1 + ^r)H^) + A(r)(l + ar)HE + je<^M^) (85) 
—jEip{-x) + je$ + jmpi/j{-x) = 0. 

Next, we need to separate out the angular dependence. This can be done by 
using the spherical monogenics [20] as in the de Sitter case. As before, we 
get the coupled equations (in matrix form) : 



with 



and 



y + C{r) jiE-e<^ + mp) + Cir) 

B = I I (88) 

j(^-e$-mp) + C(r) -y + Ci^) 

where A(r) is given by equation (80), and 
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C(r) 



A(r) 



M 



+ j{E -e^)X{r). 



(89) 



r 

In (86) Ui and U2 are the reduced functions defined as in the Schwarzschild 
case : 



Ui = ru, U2 = jrv, 



(90) 



and the primes denote differentiation with respect to r. Equation (86) can 
also be written as 



r r'^ \ U2 J 



B 



Ui 
U2 



where 



A 



1 - A(r) A(r) \ 
A(r) 1 - A(r) j ' 



(91) 



(92) 



We now look for power series solutions with a definite radius of conver- 
gence around the horizon, by introducing the series 



Ml 



k=0 k=0 

where there are now two values of 1], which are 



(93) 



(94) 



with 



r± = M ± - g2. (95) 

In the Reissner- Nordstrom case there are two event horizons, the outer (r_|_) 
and the inner (r_), not just one horizon as in the Schwarzschild case. 

We first consider 77+. As in the Schwarzschild case jH], the index s 
controls the radial dependence of if) at the horizon, which means that it 
represents a physical quantity. To find the values of s, we substitute (93) 
into (91) and set r]+ = 0. This gives 



2JM^ - g2 



r^AB 



ao 



(96) 
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where we substitute 



r+ = M + y^M2 - g2 (97) 

for r in equation (96). By rewriting equation (96) in terms of a single matrix, 
setting rj^ = 0, and setting the resulting determinant to zero, the indicial 
equation is found to be 



det 



A'B -{2 J M^-q^)-i 



= 0, 



(98) 



from which we find the solutions for the two indicial roots : 

, 2 



s = 0,--+j{E-e<^) 



(99) 



By repeating these steps for ?7_, the indicial roots which are now 



(M - - q- 



(100) 



We now derive the Hawking temperature Th- We shall need the non- 
covariant conserved current (which generates the streamlines whose curves 
are timelike), which can be expressed as j4lij : 



^ = M^7oV') det(/i)-^ 

This satisfies the flatspace conservation equation V- J = 0. 
Consider first the case with 77+ (and r+). We express t]^ as 



(101) 



Tj = exp 



+ j(^-e<f) 



M + - q' 

VM^ - q^ 



ln(r-(M+^M2-g2)) 
(102) 



We are now able to write 



ln(r - (M + ^/M'^ - q^)) = In | (r - (M + ^M'' - q^)) \ 

+ j arg(r - (M + y/M'^ - q^)) 



(103) 
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with the choice of argument 

arg(r-(M + yM^)) = | r_'^< r < r+. ^^^^^ 

Before proceeding further, we spht E into real and imaginary parts as 

E^Er- je. (105) 

If we now take the hmit r — > (M + -\/M^^^^) from above and below, we 
find that the 70 component of J is given by 



70- J = Bi{9, 0)e-2^* I (r - (M + ^M'' - q')) 
1, 



, ,, ' e) 



exp 



27r 



\ ; 



{Er-e^) 



r > r+. 



, r_ < r < r+. 



(106) 



where Bi{6,(f)) is a positive-definite finite term. 
Next, we calculate the radial component of J as 



( MJrJ M'^-tp-\ 



e'^-J = S2(^,(/))e-2^* I (r- (M + x/M^^^)) P V^?^ 

'1, r > r+. 



X 



— exp 



27r- 



-{Er-e^) 



(107) 



, r_ < r < r+. 



where B2{0,(f)) is a positive-definite finite term. 

From equation (107) we arc now able to derive the Hawking temperature. 
By taking the ratio of the outward flux to the total flux. 



1 



^^iM±^^^Er - e$) 



we identify a Fermi-Dirac distribution with temperature 



(108) 



+ 1 
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H 



2nk 



B 



(109) 



As a comparison, see, for example, jSI]. When g = 0, we reproduce the 
Hawking temperature as that of the Schwarzschild case [H] . 

By repeating the calculations for ?7_ (and r_), and with the following 
choice of argument (see also, for example, [IHj ) 



arg (r - (M - - g2)) 

the Hawking temperature emerges as 



0, r_ < r < r+, 
—71, r < r_, 



110) 



1 



H 



27xkf 



M^-q^ 



M + JM^ - g2 
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which is negative ! This is in fact impossible (see, for example, [SUliniEl)- It 
is actually the case as far as an observer outside the black hole is concerned. 
But, from the viewpoint of an observer inside the Cauchy (i.e., inner) horizon, 
the Hawking temperature of this horizon is positive ! In this sense, the 
Cauchy horizon behaves like the event horizon of the de Sitter universe. So 
it is also unstable. But, in this case, instead of expanding, it tends to collapse 
pij . If it were expanding, it would collide with the outer event horizon and 
the black hole could explode; black holes would then be very rare if exist at 
all. 



4.2 Klein-Gordon equation on Schwarzschild black hole 
metric 

We now discuss another application, the Klein-Gordon equation in the Schwarzschild 
black hole background (see also |SE])- We still employ the symbol j to rep- 
resent right-sided multiplication by la^ and work in the Kerr-Schild gauge. 
The Klein-Gordon equation is |SHl HI IHZl IHS] 

{V^ + ml)ij = (112) 
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with V = h{e'^)dn + e^a;(e^) x, and we use the separabihty ansatz 

^P = u{r)Q{d)e^""''e~^^\ (113) 

By substituting %l) into equation (112) and separating out the equations for 
e(^), we get 



u + 



-AjEMr + 2r - 2M 
r(r - 2M) 



Co 



M = 0, 



(114) 



where the primes denote differentiation with respect to r, and Ci and C2 
are separation constants. In this case we have 2A'/Ci — —2ME^ + C2 = 
—2ME + jE + where the further constant C appears in the spherical 
harmonic equation for 0, 



smO ' ^ sin^ 



)e = o, 



(115) 



by which we get C = /(/ + 1). Note that we have used A(r) = M/r here. 
Again, we introduce the series 



k=0 



(116) 



to find the power-series solutions around the horizon, where 7] = r — 2M. By 
putting this series into equation (114), we find the two indicial roots 

s = 0,AjME. (117) 
The non-co variant conserved current in the Klein- Gordon case (with h = 1) 

h(lm(-ijj*h(V) 



IS 



J 



det/i"^ 



rrir. 



(118) 



Note that deth = 1 in the Kerr-Schild gauge. We will discuss the case for 
s — AjME (s = is trivial). We first write r]^ as 



T]' = exp[(4jME) ln(r - 2M)]. 



(119) 



Write now 



ln(r - 2M) = In | r - 2M | + j arg(r - 2M) 



(120) 
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with the choice of argument 

As before, we spht E into real and imaginary parts as E = — je. 

By taking the hmit r — > 2M from above and below we find that the 70 
component of J is given by 



,„.J ^ M^Je-"-'im) I r-2M |-»'. x { ^^^^j^^^^^j^ ^ > 

(122) 

where A{9, 0) is a positive-definite finite term. The radial component of J is 



e'--J = S(e,0)(^)e-2^*|^-2M|«^^x| ^' , ^ > (123) 

^ rup' ' ' \ exp[87rM£;^], r < 2M. ^ ' 

where B(Q^ 0) is a positive-definite finite term. Upon taking the ratio of the 
outward flux to the total flux 



(124) 



we identify a Bose-Einstein distribution with Hawking temperature 



4.3 Klein-Gordon equation on Reissner-Nordstrom black 
hole metric 

We now study the Klein-Gordon equation with a Reissner-Nordstrom back- 
ground. Following the same steps, we find 



-4j{E - e$)Mr + 2r - 2M + 2j{E - e<^)q^' ^/ ^ ^ q (^36) 



^2 _ 2Mr + g2 



where C is a constant and defined via the spherical harmonic equation as 
above. 
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Introducing series as in equation (116), we get two values of 77 as before. 
We first consider rj^ (for r_|_ = M + ^/M'^ + q^). By substituting the series 
into equation (126), we find the two indicial roots 



s = 0,j(^-e$) 



(127) 



From equation (118) together with equations (103), (104), and (105), by 
taking the hmit r — > (M + \/M'^ — q^) from above and below we find that 
the 7o component of J is given by 



J = B^{e, m—)e-''' I (r - (M + - q^)) 



7o 



,(-1+2- 



f+\/Af2-q^)^ 
~7W^ X 



f 1, 



exp 



27r 



r > r+, 



{Er - e$) 



, r_ < r < r+. 



where Bi{9, 0) is a positive-definite finite term. 
The radial component of J is 



;i28) 



e^-J ^B.2{e,(t>){—)e 
nip 

r 1, 



^ ^^-2e^ 



(m+\/-M2-92^'^ 



exp 



27r 



-{Er-e^) 



, r_ < r < r+. 



(129) 

where -62(^5 0) is a positive-definite finite term. Upon taking the ratio of the 
outward fiux to the total fiux. 



1 



exp 



27r 



A/M2-g2 



(£;^ - e$) 



(130) 



- 1 



we identify a Bose-Einstein distribution with Hawking temperature 
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Th 



1 



- g2 



(131) 



2nk 



( 



M + y^M2 - g2 



2 



By performing similar calculations for rj^ (for r_ = M — v^M^ — q'^) with 
the choice of argument of (110), we find the Hawking temperature 



which is negative as in equation (111). Further discussion can be found in 
0. 

4.4 Magnetic Monopoles 

In this paper we add magnetic monopoles to the discussion of the previous 
section to generalise the corresponding expression for the Hawking temper- 
ature jlHll2ll2nill2lllSlllllE3EZ|- The magnetic monopole, which is a 
point magnetic charge, was introduced by Dirac to remove the asymmetry 
in Maxwell's electromagnetic equations. Dirac derived a quantization con- 
dition for the magnetic charge : if magnetic charge exists, then it will take 
discrete values. In Dirac's proposal, the existence of any magnetic monopole 
is responsible for the the quantization of electric charge. With h = c = 1, 
the Dirac quantization condition for electric charge e and magnetic charge 
Qm (such that B = qu/'^'^'^r for a point charge) is 



where N is an integer. 

We can express the vector potential A for a magnetic monopole by the 
following equation : 



The 'string' along ^ = is unobservable provided that is an integer. 
The additional term due to monopoles causes the change + q\j (see 

[17j for details), so the A(r) changes slightly to 




(132) 



2eqM = N, 



(133) 




(134) 
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r 2r 

All other terms stay the same. We now make a different ansatz 



Kr) ^ ^ - (135) 



ijj = u{r)ijj^ + Iv{r)ij~ (136) 
where ip'^ and ip" are given by 



yvnk 



V;+ = ^ + 1 + + V/ + 1 - A:) +$[ 

+(V^ + 1 + A; - + 1 - A;) ' 

and 



(137) 



-(V/ + 1 + A; - V/ + 1 - A;) ' 



(138) 



with eigenvalues ((/ + 1)^ — k"^) ± ^J{1 + 1)^ — A;^ = ± 1) respectively, 
where +<I>J"*^ and ~<l>^*^ are 



-^J"'^ (X ±^Z(Z + 1) - m{m - l)Yi"''' - Ia2{l - m)Y/^+^\ (139) 
and 



-^r'' ^ ±(^ - m)r,'"'= + Ia2^l{l + 1) - m(m - 1)1^"^+^ ^ (140) 

Details of the above monogenics are given in Appendix A. After some algebra, 
we get the coupled equations as usual as 



with 



and 
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B 



where 



/ ^ + CW j{E-e^ + mp) + C{r) 

[j{E-e^-mp) + C{r) + CW 



X(r) M 



(143) 



and Ui and M2 are the reduced functions defined as in the Schwarzschild case: 



Ui = ru, U2 = jrv. 
Here the primes denote the differentiation to r, so that 



U2 



where 



,^/l-A(r) A(r) 
[ A(r) 1 - A(r) 

We look for power series solutions of the form 



k=0 k=0 

where we have two values of 1], which are 



r]±=r - r±, 



(145) 



(146) 



(147) 



(148) 



(149) 



with 



r± = M±JM^-{q^ + ql). 



(150) 

We now consider 77+. As in the Schwarzschild case jH], the index s 
controls the radial dependence of ip at the horizon, so that it represents a 
physical quantity. After some algebra we get the following indicial equation 



det 



A'B -{2 J M^-{q^ + ql))^J 



0, 



(151) 
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from which the solutions for the two indicial roots are 



s = 0,--+j(^-e$) 



(M + ^M2-(g2 + g2^) 



M2 



(152) 



In deriving the Hawking temperature Th we shall need the non-covariant 
conserved current (which generates the streamlines whose curves are timelike) 



(153) 



J = 4(^70 V^) det{h)-^ 

which satisfies the flatspace conservation equation V- J = 0. 
For the solution 77+ (with r+), we can express r/"* as 



/ 



7] = exp 



We can now write 



-i+j(^-e<f) 



X ln(r - (M + ^M^-{q^ + ql,)))} . 



(154) 



ln(r - (M + JM^ - (g2 + ^2^))) = In | (r - (M + - (g2 + 4^))) 



with the choice of argument 



+ j arg(r - (M + y'M^ - (g^ + g2^))) 

(155) 



arg (r - (M + ^M^ - {q^ + ql,))) = | 

We split E into real and imaginary parts as 

E = Er- je. 



0, r > r+ 

— 7r, r_ < r < r+. 



(156) 



(157) 



If we now take the limit r — >■ (M + yM^ — (g2 -|- g|^)) from above and 
below we find that the 70 component of J is given by 
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70- J = B,{e, 0)e-2^* I (r - (M + ^ - {q^ + ql))) 

1, 

2 



X 



exp 



(E, - e$) 



where Bi{6,(j)) is a positive-definite finite term. 
The radial component of J is 



r > r+, 
, r_ < r < r+, 
(158) 



)e-2rf I 



f 1, 



(M+ jM2-(g2 + ^2^))) I 



X < 



exp 



27r^ ^= ==^(Er - 6$ 



r > r+, 
, r_ < r < r+, 

(159) 

where B2{9, (p) is a positive-definite finite term. 

We arc now able to derive the Hawking temperature. By taking the ratio 
of the outward flux to the total flux, 



+ 



exp 



27r 



(M+^/m2^)2 



{Er - e$) 



(160) 



+ 1 



we identify a Fermi-Dirac distribution with the revised Hawking temperature 
(for r+) 



- (g2 + ql) 



M + ^M^- (g2 + ql) 



(161) 



4.5 Electromagnetism in the Schwarzschild black hole 
background 

We first find a covariant Maxwell equation in the presence of a gravitational 
field. We extend the 'wedge' equation 
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V Ah{a) = kS ■ h{a), (162) 

to include higher-grade terms jH]. To keep the derivation general we include 
the spin term, so that 

VAt{aAb) = [V At{a)]Ah{b) -h{a) A[V At{b)] 

= 4h{a) ■ S] A h{h) - Kh{a) A [h{b) ■ S] (163) 

= kS X h{a A b). 
More generally, we can write 

V A h{Ar) = h{V A Ar) + K{Sh{Ar))r+l. (164) 
We can therefore replace V A F = by 

VAJ^-kSxJ^ = 0. (165) 
Now we use the rearrangement 

V ■ (/i(J^) det(/i)-i) = IV A{Ih{T)det{h)-^) 

= IVA{r\lJ^)) (166) 

= ihr\v A{ij^) + X s], 

to write 

V ■ [M(F) det(/i)~^] = J, (167) 

with F = V A A, where A is the vector potential and J is the charge current, 
as 

V-J^-kS-J^ = Ih{JI) = J. (168) 
By combining the equations above, we get 

VJ" - KSr = J. (169) 
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To find the free-field stress-energy tensor, we apply the definition 

det(A)^(„)(^mdet(A)-') = Tr'(a), (170) 

to get 

Temh-'ia) = |det(/i)%(„)(MF)MF)det(/i)-i) 

(171) 

= h{a- F)- T -\h-\a)T ■ T. 

Hence, 

%rn{a)^^h{a)■F)■J'--aT■T ^{a-T)-J'--aT-T ^--TaT, (172) 

which is the natural covariant extension of the gravitation-free form —^FaF. 
This tensor is symmetric, as expected for fields with vanishing spin density. 

We are now ready to investigate electromagnetism in a Schwarzschild 
black hole background. Here we use the Kerr-Schild gauge 

^(a) = A(r)a-e_e_, (173) 

where 

, . M , , 

A r = — , 174 
r 

and e_ = 70 — e,.. Wc begin by writing the electromagnetic field J- (see 
Appendix B for details of the vector spherical harmonics) as 

jr ^ F{r, i)|(l + (7,)Xr + G{r, 1)^(1- (7,)Xr + H{r, t)arY{^ (175) 
where the functions F, G, and H are complex, and 

V(a,Xr) 

v(a,rr) 



1 \Ki + 1) 

= ^ rYr^ (176) 

1 Jl{l + 1) 

= -X™+-^^^ ^f^rl^r, (177) 

r r 

1 JZ(Z + 1) 
= -X^+-^^^- -Yl^. (178) 
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Here we have defined the normahsed harmonic 

1 



yi(7Ti) 



where 



(179) 



L = -X X V. 



We also know that 



2/ -I- 1 (] — rnV 
with Yf"^ = {-1)'^{Y{^)* . By forming VT = 0, we get 



(180) 



(181) 



/ 


l+2f 








\ 







-1 















1 




V 











/ 



G I + 



I 



1 

r 



2M o\ 

-1 

\) 

r 

v^/(/ + i) 



2M\ VHi±ll 



/(/ + 1) 
5? 



V 



I 2M^ V^(^ + i) V^(^ + l) 2 
V r ) 2r 2r r 



G 
H 



J 



(182) 



Upon separating of variables as F(r)e G{r)e and H{r)e the 
radial part of the field emerges in terms of the transverse components as 



H 



2jrE 

The radial equations then become 



1 - 



2M 



r / 



)f + g) 



(183) 
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/ 2M\ , F f 2M\ 1(1 + 1)// 2M\ \ 

(184) 

Next, series expansion of these equations, as before at the horizon using 
ri = r — 2M, yields 




TEW 

V 

The roots of the indicial equation are 



s = 0, -1 + ^jEM. 
The regular solution at the horizon satisfies 



G 



(186) 



= 



(187) 



(188) 



r=2M 



As before, we now proceed to calculate the Hawking temperature. Since the 
imaginary part of the non-zero root in this case is the same as in the Dirac and 
Klein-Gordon cases, we expect to recover the same Hawking temperature. 

If we look at the absolute values of the real part of the non-zero roots in 
the Klein-Gordon, Dirac, and electromagnetic cases discussed in former chap- 
ters, we find that they are equal to 0, |, and 1, respectively. These are also 
the values of the spin of the particles involved. Moreover electromagnetism, 
as discussed here, is completely classical. It seems that geometric algebra 
can detect that Maxwell's theory concerns spin-1 particles (i.e., photons) 
even though the picture of bosonic particles carrying forces is a quantum 
description! The quantal 'geometric' nature inherent in the theory can in 
some sense be recognised by geometric algebra! 

We now calculate the Hawking temperature. We use the stress-energy 
tensor 



35 



Te„(a) = -]^Ta:F. (189) 

We first write rf as 

rf = exp[(-l + AjEM) ln(r - 2M)], (190) 

and then write 

ln(r - 2M) = In | r - 2M | + j arg(r - 2M) (191) 
with the choice of argument 

As before, we spht E into real and imaginary parts a.s E = Ej. — je. Then 
we calculate 7^m(e'"), which reads 



TUe') = Z(0, ^)e- \r-2M x { ^i^p.^.^^j^ : J (193) 

where Z(6', 0) is a finite term. By taking the ratio of the outward stress- 
energy tensor to its total : 



em+ 



T — T 

-'em+ ■'■em- 



we identify a Bose-Einstein distribution with Hawking temperature 



■H 



(194) 



(195) 



5 Axially Symmetric Black Holes 

5.1 Klein-Gordon equation on Kerr black hole metric 

We now begin the calculations using GTG [HI HH]. We shall discuss the 
Klein-Gordon equation in the Kerr black hole background As usual, we 
will employ the symbol j to represent right-sided multiplication by la^, and 
work in the Kerr-Schild gauge. The Klein-Gordon equation is 
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(P^ + ml)ij = (196) 

with V = h{e'^)d^ + e^uj{e^) x, and using the following Kerr-Schild gauge 
[TH] (as a comparison, see also [121 El) : 



hia) = a + di{r, 9)a ■ e_e_ + (tt — l)a. A (6^64)6^64 — ■ ete^ 

+ -7^ sin 9a ■ e+e- + -^a ' ' 



where 



and 



with 



Mr 

^i(r,^) = — (198) 

p2 



p2 = + cos^ ^, (199) 



e_ = 7o - Cr, e+ = 7o + e^. (200) 
We now employ the separation ansatz 

^ = u(r)e(^)e^'"^e-^'^*. (201) 

By substituting into the Klein-Gordon equation and separating the equa- 
tions for 0(6*), we get 



C 

{r'^-2Mr+L'^)u" + {2r-2M+2jmL-4jEMr)u+{Fr{r)-—)u = 0, (202) 



where C is a constant, m is the azimuthal quantum number, the primes 
denote differentiation with respect to r, and Fr{r) is given by 

Fr{r) = -2] EM + 2rME'^ + r'^{E'^ - ml). (203) 
The equation for 0(^^) is given by 

sin^ OQ" -2 cos OQ' + {Fe{e) + C)e = 0, (204) 
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where = 



d 

ab- 



solutions are spheroidal harmonics in the standard form 



SmniLp, cos6) with p = E — rrip, and C is the same constant as in the 
previous equation (C = \mn) (see |lj for details), and Fg{9) is given by 



Fe{e) = {E' -m;)L' cos' e 



m 



(205) 



In this case we need the equation for u{r). Again, we introduce the series 



u 



k=0 



(206) 



to study the solutions around the horizon. We find two values of f] as in the 
Dirac case. First we consider the case of r]_^_ (for r_|_ = M + y/M'^ — L^). By 
substituting the series into the equation, (202), we obtain the two indicial 
roots 



0,J 



E ({^M + VM^-L^y + lA - ml 



v/M2 - L2 



(207) 



From we know that Qh = L/{r']_ + L^), where VLh is the angular velocity 
of the horizon of the black hole. So we can rewrite the two indicial roots 

above as 



s = Q,j{E~mnH) 



M + y/M^ - m +L 



(208) 



The non-covariant conserved current in the Klein-Gordon case (with h 
1) is given by 



J 



det h 



Consider first the case for rj^ (with r+). We express r/"* as 



exp <^ \3{E - rnQn) 



ln(r - (M + - 



We can now write 



(209) 



(210) 
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ln(r - (M + - L^)) = In | (r - (M + ^f^P^^)) \ 

+ J arg(r - (M + VM^ - L^)) 

with the choice of argument 



(211) 



arg (r - (M + ^ - L^)) = | 



0, r > r+ 

< r < r+. 

Before proceeding further, we spht E into real and imaginary parts as 



(212) 



E = Er- je. 



(213) 



If we now take the hmit r (M + V — L^) from above and below, we 
find that the 70 component of J is given by 



E.. 



7o-J = i^i(^,0)(— )e-^^*x 



(r- (M + VM23X2)) 
1, 



('m+Vm2-l2)%z,2 ' 
,(-1+2^^ ^^^^^ ^) 



exp 



27r 



r > r+, 

(M+VM2-L2)^+L2)(E^-mnjf)' 

■ ■ , r_ < r < r+. 



VM2-L2 



where Ki{9,(j)) is a positive-definite finite term. 
The radial component of J is 



(214) 



E, 



e^-J = ir2(^,0)(— )e-2^*x 



M+\/ m2-l2)^ + l2^ 

r - (M + - L2)) 1^^^^ VM2-i2 ) ^ 

f 1, 



exp 



27r 



VM2-L2 



r > r+, 
, r_ < r < r+. 



where ^^2(^,0) is a positive-definite finite term. 

Upon taking the ratio of the outward flux to the total flux. 



(215) 
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1 



exp 



27r 



(216) 



we identify a Bose-Einstein distribution with Hawking temperature 



T, 



H 



2TrkB 



(217) 



The same results can also be found in jHO], jHHHl, jHE]; [Hj; and [TT] . 



By running through the calculations for ?]_ (for r_ = M — \/M2 — L^), 
and taking the choice of argument as in the Dirac case, with indicial roots 



s = 0, -J 



(m + -L2)^ + L^") (E - m^]^^ 



VM2 - L2 



(218) 



we find the Hawking temperature is 



27rA;f 



\/M2 - L2 



M + VM2 - L2 



(219) 



We can also do the analysis for the Kerr solution in oblate spheroidal 
coordinates JF;. We take the gauge fields in Kerr-Schild form ^Hj : 



h = a + a -11, 

where / = \/a{l + 11)70. In oblate spheroidal coordinates 



(220) 



{x, z} = {L cosh U cos V cos 0, L cosh U cos V sin 0, L sinh U sin V^} 
we have 

1 



n 



L cosh U 



{eulo + L cos Va^ 



(221) 
(222) 



40 



and 



MsinhC/ 



(223) 



L(cosh2 U - cos^ V) ■ 
The horizons of this solution are given by 



Lsinh [/ = M ± - 



(224) 



With the separation ansatz 



ip = U{U)V{V)e 



,jm<p -jEt 



(225) 



we successfully obtain separated equations for V{V) and U{U). These are 



The roots of the indicial equations are the same as those above. Note that 
r — L sinh U. 

5.2 Dirac equation on Kerr black hole metric 

Wc do not discuss Dirac equation on Kerr black hole metric in this paper. 
We only noted that the roots of the indicial equation obtained for that case 
are given by 




(226) 



where C is a constant, jjl = sinT^, and 




(227) 





(228) 



- L2 
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where m' = (m + |) and VIh = L/{r\ + L^). This gives the same Hawking 
temperature as in the Klein-Gordon case. The temperature of the black 
hole is determined by the imaginary part of the non-zero root. In this case, 
the root is the same as in the Klein- Gordon case apart from the difference 
between m' and m, which does not affect the temperature. 

We do not discuss the Kerr- Newman case in this paper since it is straight- 

forward. The only change is to send M ^ M — 

6 Conclusion 

The results in this paper have shown us that geometric algebra, and specifi- 
cally spacetime algebra, substantially simplifies the study of the Dirac equa- 
tion and its consequences. It provides a powerful and flexible tool even for 
the analysis of strong gravitational backgrounds for which the Gauge Theory 
of Gravity (GTG) was also employed. 

We have studied how geometric algebra and GTG can be applied to 
black holes. We have considered the Klein-Gordon (spin-0), Dirac (spin- 
1/2), electromagnetic field (spin-1) equations in the Schwarzschild, Reissner- 
Nordstrom, and Kerr black hole backgrounds. The Hawking temperatures 
were also derived for all cases and compared with the results of 'standard' 
calculations. Remarkably, they agreed. We draw the conclusion that we can 
use the first-quantised quantum theory instead of the second-quantised one 
originally used by Hawking [30 in treating the black hole evaporation 

In this work it was found that the absolute values of the real part of the 
non-zero roots in the Klein-Gordon, Dirac, and electromagnetic field cases are 
the values of the spin of the corresponding particles yet electromagnetism, 
at least, is completely classical. It seems that geometric algebra can tell 
that Maxwell's purely classical theory concerns spin-1 particles (photons) 
even though the picture of bosonic particles carrying forces is a quantum 
description! In some sense, the quantal 'geometric' nature of the theory is 
recognised by geometric algebra! 

The result for the 'Hawking temperature' of the de Sitter cosmological 
horizon is also interesting. If a thermodynamic system has a negative tem- 
perature then that system should be unstable, and tend to expand or col- 
lapse. But, if we instead abandon the concept of particles as being observer- 
independent in this case, we find the usual positive 'Hawking temperature'. 
Geometric algebra was also successfully applied to treat the Unruh effect. 
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In the separation E — Ej. — jt (or E = Er + je), the e term is never 
equal to zero unless the potential of the event horizon is infinite. In the cases 
discussed here the potential is finite, so that e is non-zero. There is not then 
any problem with the density, which is not normalizable at the event horizon 
if e is equal to zero; it is normalizable even at the event horizon. Moreover, 
even if -0 is singular at the event horizon, the normalization integral is finite 
because the density is normalizable there. We may conclude that there is no 
problem with the normalizability. 

We compare this situation with the tunnelling effect through a finite 
potential barrier in quantum mechanics. The particles in a box defined by 
such a potential will leak out from the box. This means that the imaginary 
component of E is non-zero, because the fiux of particles tunnelling through 
the barrier is non-zero. Recall that the event horizon of a black hole is the 
boundary which the particles which come from inside the black hole cannot 
pass through unless they have superluminal velocity. But this does not mean 
that the potential of the event horizon itself is infinite. The potential remains 
finite, and from the viewpoint of quantum mechanics (which also supports 
the idea of non-locality) it can be penetrated. That is why e should be 
non-zero in this case. 

According to second-quantised theory, the source of particles in Hawk- 
ing radiation comes from somewhere very near to, but outside, the event 
horizon of a black hole, and is due to vacuum polarization. Particles with 
positive energies escape to infinity and those with negative energies fall into 
the black hole. The black hole therefore 'evaporates'. This argument cir- 
cumvents superluminal effects (the breakdown of causality) at the expense 
of introducing the vacuum polarization near the event horizon. If we accept 
the non-locality inherent in Quantum Mechanics and the possibility that 
causality breaks down in very strong gravitational fields (in which we should 
employ a quantum gravity approach), we can use the first-quantised theory 
instead of second-quantised theory when dealing with problems involving 
sources with very strong gravitational fields, such as black holes, as in the 
evaporation process studied here. 

A hard problem in black hole physics concerns the endpoint of black hole 
evaporation. The difficulty is that at this late stage the effect of quantum 
gravity can no longer be ignored. This is a future topic for research within 
the GTG framework, which has proved very powerful in this work. It will be 
interesting to see what the application of GTG will say about black holes as 
well as quantum gravity in the future. 
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7 Appendix A 



Assuming that spherical monogenic ipj^ is an eigenstate of the x A V and J3 
operators (see [20] for details), we find after some calculation that 



- X A = ^ > (229) 

J3^[" = (m + i)V^[" -1 - / < m < / (230) 

The highest-m case satisfies 

J+^l = 0, (231) 
and, after introducing a convenient factor, this is solved by 

ipl = {21 + l)P/(cos0)e''^^"^ (232) 
where the associated Legendre polynomials P[^{x) are given by 

which satisfy the following recursion relations : 

(1 _ x^) \ ^ ' + mxPr(x) = -(1 - x^f^Pr+\x), (234) 
ax 



{1-x Y ' -mxPrix) = {l-x^ f\l + m){l-m + l)Pr'\x). (235) 
The lowering operator J_ acting on produces 



J_V = [~deiJ + cot ed^tljla3]e-'^^''' - ^^(^ + as^ag), (236) 
and the derivatives applying on ipl give 

[-deipl + cot ed^i)\la^]e-^^''^ 

= {21 + l)[-dePlicose) - /cot^P/(cos^)]e('-i)'^^"3 (237) 

= {21 + l)2/P/-^(cos^)e('^l)'^^"^ 
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Since 



(T^ = (J2e'^'''^ (238) 
we get after some manipulation the following spherical monogenics : 

^™ =[{l + rn + l)Pl^{cos e) - P;"+\cos 9)Ia^]e""^^''' , (239) 

where / is a non-negative integer, m runs from — (/ + 1) to / and P[\x) are 
taken to be zero if \m\ > I. 

In order to get generalise the spherical monogenic functions suitably, we 
need to generalise the spherical harmonics by including dependence on the 
third Euler angle x- We begin by considering the quantum treatment of 
rotating bodies, which has an obvious relation to the generalised spherical 
monogenics. We express the space axis angular momentum operators in 
terms of the usual Euler angles ^, 0, x : 

jL^/h = — sin (pdg — cot 9 cos (t)d^ + coscc^^ cos 09^, 

jLy/h = cos 4>dQ — cot 9 sin 09^ + cosec6' sin (pdy., (240) 

The corresponding form for the components with respect to the body axis is 
given by 



jLyh = sin xde — cosec^ cos xd(j, + cot 9 cos xd^ 



X' 



jL\/fi = cos xd$ + cosec9 sin xO(j) — cot 9 sin xO^, (241) 
These operators satisfy the commutation relations 



[LI, LI] = jfiLl- [L\,L\] = -jhL\- [LI, L\] = 0, etc., (242) 

so that the space and body operators commute with each other, but their 
commutators amongst themselves differ in sign. 

In geometric algebra, we need a coordinate-free translation of the angular 
momentum operators in terms of the generators of an active rotation, which 
will be described by a rotor R [RR = 1). The effect of R on any multivector 
M is 

M ^ RMR. (243) 
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The configuration of a rotating body is described by a rotor R giving tfie 
directions of the body axes in terms of the space axes 0"^ : 

Cfc = RakR. (244) 

Let us imagine a molecule as our rotating body. We can set up a wavefunction 
ip{R) on this configuration space. Note that, under an infinitesimal space 
rotation in the plane of bivector B, 

ilj{R) ^ i/j{e'^/^R) « V(^) + -{BRdR)^; (245) 

where Br is the multivector derivative with respect to R. The 'space' com- 
ponent of the angular momentum operator for rotations in the plane of the 
bivector B is therefore 

L% = ^jhiBRdn). (246) 

It is easily verified that 



[L%L%] = C2j^f{{ARdR){BRdR)-{BRdR){ARdR)) 

= C^jhyHARdnB) - {BROrA)) (247) 
= -Ih^liBA - AB)RdR) = -jhLjy^B. 

The sign of these commutators reflects the fact that we have chosen to work 
with bivector components rather than the usual vector components. 

The 'body' components are fixed in the rotating frame, so we send B 
RBR in the previous formula, and define 

4 = ^jh{RBdR). (248) 
These operators commute with the space operators : 



[L%L'b] = i^,jh)^{{ARdR){RBdR)-{RBdR){ARdR)) . . 

= {yhy{{ARBdR)-{RBdRA)) = 0. ^ ' 

A similar calculation yields 

(250) 
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as expected. 

We are now ready to resume seeking the generalised spherical monogenic 
functions. The unnormalised generalised spherical harmonics related to the 
functions we need will be written 

Yr\li. 0, x) oc Pr\^^) exp / W + kx). (251) 

where Pl^^{^) are the generalised associated Legendre functions. The Y";"**^ 
are eigenf unctions of the space axis and body axis angular momentum oper- 
ators given by 

LlYi""^ = mYl^^ , Llrr^ = kYi""^ , LlLlYi"'^ = l{l + l)Yl^^ . (252) 
The Clifford versions of these operators are 

Ll^-^-{Ia^RdR), Ll^-l{RIcr,dn), (253) 

where j is the Clifford version of the bivector imaginary operator. The solu- 
tion is standard and P[^''{^) is given by pOj 

pr'if^) = ^(1 - f,p/\i + fir^/'pjT'^Kfi), (254) 

where a = \m — k\, P = \m + k\, n = I — a/2 — P/2, and -P^"'^-* are the Jacobi 
polynomials 



( -\\n -I m 

The allowed values of I, m, k are 

/ = 0,^,1,|...; -l<m,k<l, (256) 

which reflects the fact that a, /?, n must be non-negative integers. The 
half-integral angular momentum states have appeared already; the simplest 
is 

Yp = cos(0/2) exp(/(0 + x)/2). (257) 
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The normalising constant in this case is 

/n!(a + (3 + n)!(2n + a + /? + 1)^^^ 
2"+^+i(n + a)!(n + /3)! ) 
and the action of the ladder operators is then 



(258) 



^m—l k 



ls^Y{^k ^ ^^i^i ^ ;l) - m(m T l)Yl 

Ih^rak ^ ±^l{l + l)-k{kTi)Yr''~^ (259) 

where the overall sign depends on ni and k. 

The first set of generalised monogenic functions '^^f^^ we need are the 
eigenfunctions of J|J|, J|, and L^, with eigenvalues (/ + |)(/ + §)■ and k, 
respectively. These are the J = / + | states; we will find j = I — \ afterwards. 
Since Yl^ is already a generalised monogenic, we develop those for m < I 
using the J_ operator, and find for m < / that 



oc ±^l{l + 1) - m(m - l)Yi'^'' - Ia2{l - m)l^"^+'*^ (260) 

where the + sign is taken if m < A;. The normalisation is obvious (note that 
Ia2e^"'^'^ — lacf). The range of I, m, k is 

1 3 

/ = 0, -l-l<m<l] -l<k<l, (261) 

so that the total number of states for a given j — I -\- \ is {2j + l){2k + 1). 
The simplest example of a half-integral state is 

+$;/^'/' = cos(^/2) cxp(/a3(0 + x)/2). (262) 

Conventional Quantum Mechanics leads us to expect another set of j = 
I — I monogenics provided that I > 0. The other set of monogenics is 



-^r'' « ±(^ - m)Yl^'' + Ia2^l{l + 1) - m(m - l)Y;'"+'^ (263) 
The range of Z, m. A; is 

1 3 

/ = -,!,-,...; -/<m</-l; -l<k<l, (264) 
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so that the total number of states for a given j — l — ^is still (2j + 1) (2k + 1) . 
The simplest example of a half-integral state is the j — doublet I — ^, 



= sin(e/2) exp(7(73(-0 + x)/2) + M cos(^/2) exp(7(73(0 + x)/2). 

(265) 

From the generalised spherical harmonics we can define the monogenics 
±$^*:^ which are the eigenf unctions of the space angular momentum operator 
J* with eigenvalues j — l:k\. For the case A; = 0, there is a simple relation 

= ±Iar +*r^(^3 (266) 

between the normalised monogenics. The sign in this relation depends on 
the sign of (m + |). 

The monogenics developed above can handle problems involving the mag- 
netic moment, which involve strong spin-orbit coupling. The conventional 
wisdom of adding L -|- S = J should then be used. We therefore need to 
find suitable combinations of ^$ states that are eigenstates of the (p — eAY 
operator. Since j, m, k are still good quantum numbers, we expect to use a 
combination of and to treat the magnetic monopole case using 

the form of ijj'^ and given by 



and 



^- = ^ ((V/ + 1 + A; + Vl + l-k) +<l>r'= 
-{Vl + l + k - Vl + l-k) -<^>r+i) 



(267) 



(268) 



which have eigenvalues {{l + iy — k'^)±y {I + 1)^ — k"^ = k,{k,±1) respectively. 

Care is necessary with the notations employed because it can be difficult 
to avoid overlapping use of the same symbols. However, it is always clear 
from the context what is meant. 
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8 Appendix B 



We construct the vector spherical harmonics from the scalar spherical har- 
monics given by 



Yr{9, <p) ^ ^ ^|q:^^^(cos^)e^-^ (269) 

with Yf'^ — (— l)"*(y^"^)*. The vector harmonics arise from the action of the 
quantum angular momentum operator 

L = -JxxV = -X X V, (270) 

where x denotes the Gibbs vector product and x is the STA commutator 
product. We define the normalised harmonic 

Xr(^, 0) = J-— LYr{e, <j>), (271) 

whose 'Cliffordised' versions of the interrelations are given by 



1 \Ki + 1) 

VXf^ = -cT^XJ" - -arY/^, (272) 



1 \Kl + 1) 
V(a,Xn = _X[" + -^^-- -arYi"", (273) 



1 \Kl + 1) 
V(a^F™) = -X™ + -^^-- -Yi"". (274) 

The other harmonics we need are related to Xp, the full set being 

XT', a^xXT' (transverse), a^Y;'" (longitudinal). (275) 
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